Negative and positive hysteresis in double-cavity optical bistability in three-level atom 
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We present novel hysteretic behaviour of a three-level ladder atomic system exhibiting double- 
cavity optical bistability in the mean-field limit. The two fields coupling the atomic system experi- 
ence feedback via two independent, unidirectional, single mode ring cavities and exhibit cooperative 
phenomena, simultaneously. The system displays a range of rich dynamical features varying from 
normal switching to self pulsing and a period-doubling route to chaos for both the fields. We focus 
our attention to a new hump like feature in the bistable curve arising purely due to cavity induced 
inversion, which eventually leads to negative hysteresis in the bistable response. This is probably 
the only all-optical bistable system that exhibits positive as well as negative bistable hysteresis in 
different input field intensity regimes. For both the fields, the switching times, the associated critical 
slowing down, the self-pulsing characteristics, and the chaotic behaviour can be controlled to a fair 
degree, moreover, all these effects occur at low input light levels. 



All-optical bistability has remained a focus of research 
for more than four decades [ll-Wand apart from its po- 
tential application as a switch [5] in optical communi- 
cation technology it has continued to remain a test bed 
for fundamental research related to cooperative phenom- 
ena @, 0] as well as to the study of nonlinear dynamical 
aspects such as self-pulsing, instabilities and chaos [1- 
1 1 II ] . A deeper understanding of issues related to quan- 
tum aspects like entanglement and cooperative behavior 
in presence or absence of instabilities would be critical 
in realizing a functional quantum computer. Due to re- 
cent developments related to cold-atoms in optical lat- 
tices [IH and atomic chips [l3|, the aspects related to 
cooperative phenomena have become vital (lij . Design 
of smaller trap size (< A) demands a treatment that al- 
lows for cooperative effects, and possibly at multiple fre- 
quencies. Cooperative effects at multiple frequencies can 
be realized in multi-level atoms, and we show the on- 
set of instabilities in such systems even at low light lev- 
els. Multilevel atoms have been recently used to create 
and control cooperative effects, such as in the multipar- 
ticle dark states [l5[ and Rydberg blockade effects [l6| in 
cold atoms in a trap. In order to understand interplay 
of multi-colored cooperative effects we explore the semi- 
classical dynamics of three-level atoms interacting with 
two fields coupling two adjacent transitions at which they 
exhibit cooperative behavior simultaneously. The atomic 
level structure itself provides the coupling between the 
two distinct cooperative branches. 

In an early work, one of the authors (HW with G.S. 
Agarwal) had shown the possibility of control of optical 
bistability [l7| in three-level atomic medium. In that 
configuration, the optical bistability exhibited by one 
(probe) field is controlled by another (control) field cou- 
pling an adjacent transition. The field exhibiting bista- 
bility experiences conventional cavity feedback, whereas 



the control field is held constant without feedback. In or- 
der for the probe field to exhibit optical bistability (OB) 
it not only needs to interact strongly with the active me- 
dia but it must experience sufficient feedback (thus sat- 
isfying the conditions necessary to exhibit cooperative 
behavior), which is obtained through an external unidi- 
rectional ring-cavity. Such configuration allows effective 
engineering of various characteristics of OB, including 
tailoring of thresholds, changing On/Off intensities of the 
output field and obtaining multistability [l8[ . These ef- 
fects were also experimentally realized [19]. Furthermore, 
three-level atom with a control field (without feedback) 
has been shown to exhibit instability in the context of 

ob m. 




FIG. 1: (Color online) (Top) Schematic of the three-level lad- 
der system interacting with the two fields Ei and E2 at fre- 
quencies wi and L02, respectively. (Bottom) The OB set up 
with the collection of atoms (within the length L) that in- 
teract with the two fields in two independent unidirectional, 
single mode ring cavities. 
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In this work, we again have two fields that couple two 
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adjacent transitions in a three-level atom, however, both 
the fields experience feedback via two independent, sin- 
gle mode, unidirectional ring cavities, as indicated in 
Fig. [TJ This system provides independent control over 
the two cooperative parameters. The fields are chosen 
to be counter-propagating within the active medium in 
order to minimize the two-photon Doppler broadening 
in the ladder-system, however our calculation is under- 
taken for a homogeneously broadened atomic gas. Such a 
configuration leads to negative hysteresis apart from the 
conventional positive bistable hysteresis. As mentioned 
above this single mode two cavity system leads to a vari- 
ety of nonlinear dynamical effects including self-pulsing 
of the two fields at their distinct frequencies and a period- 
doubling route to chaos in the lower cooperative branch. 

Nonlinear dynamics related to OB has traditionally 
been associated with inclusion of multi-modal treatment 
of the cavity field 0, [Hj], here, we have just two modes 
corresponding to the two monochromatic fields coupling 
the adjacent transitions. The theory of two-photon am- 
plifiers and absorbers forms the traditional basis of study- 
ing these systems, which traditionally reduce the problem 
to an effective two-level model and has been extensively 
studied [22|, (23|. An early work describes the possibil- 
ity of obtaining chaos for atoms interacting with a sin- 
gle mode, where chaos occurs in the upper-branch for 
extremely large cooperative parameters accompanied by 
large atomic and cavity detunings [241 ]. There have been 
other studies relating to the dispersive regime of two- 
photon OB which deals with the effects of cavity detun- 
ing that controls both the fields simultaneously (25|. In 
our model, the system exhibits chaos at sufficiently low- 
light levels in the lower branch of the OB response with- 
out being restricted to any special regime. Moreover, 
the independent cavities allow for an effective control of 
the feedback of the fields, thus tailoring the cooperative 
character. It should be noted that all the rich nonlinear 
dynamical features arise purely due to the interplay of 
the two cooperative branches. 

Apart from a plethora of nonlinear dynamical features, 
the fields exhibit cavity induced inversion and positive 
as well as negative hysteresis OB, as one varies the in- 
put intensity. To the best of our knowledge, we are not 
aware of any other all-optical bistable system that ex- 
hibits chaos at low-light levels or even the negative and 
positive hysteresis OB. Our results are also in contraven- 
tion with those reported earlier, in particular we observe 
enhanced output field resulting from creation of popula- 
tion inversion between the states |1) and |2), in presence 
of concurrent feedback for fields coupling both the tran- 
sitions. The effects are particularly significant as they 
occur at low light levels and thus could be used as a 
nonlinear component in optical circuitry. The details of 
the nonlinear dynamical studies are presented in detail in 
the companion paper [26j j . The atom-field density matrix 
equations and the field equations governing the dynamics 



of the OB system in the mean-field limit are given as 



-gj- = Kx [-xi(l +iOi) + y% + 2iCipi2,] 

-gj- = K 2 [-X 2 (l+i9 2 ) + y 2 + 2iC 2 p 23 .} (1) 

These equations describe the three-level atom coupled to 
the two fields at frequencies co\ and ui 2 , which experience 
feedback through two independent cavities characterized 
by the cavity decay n\ and k 2 , with cavity detuning 5\ 
and 8 2 (scaled as 9i = Si/ Ki where i = 1, 2) respectively. 
The total Hamiltonian in the dipole approximation af- 
ter undertaking the rotating wave approximation is given 
as: H = ft(Ai + A 2 )|1}(1| + SA 2 |2}(2| - (d 12 .E 1 \l)(2\ + 
<i 2 3-£ : 2 12) (3| +h.c), where Ai and A 2 are the atomic de- 
tunings. The detailed density matrix equations and the 
scaled variables along the lines of Ref. Q are all explicitly 
given in Ref. [26j . 

We consider the mean-field limit wherein a single pass 
through the ring cavity only marginally affects the fields, 
and the strong cooperative nature arises due to the ex- 
tremely large photon lifetime of the fields within the cav- 
ity, as the transmission coefficient of the cavity mirrors 
is chosen to be negligibly small. The last two equa- 
tions in Eq. (JXJ) arise due to the cavity feedback and 
the fields are determined self-consistently by the coupling 
with the three-level atom via the cooperative parame- 
ters (Ci and C 2 ). Under steady state conditions one can 
obtain the various domains of stability for different in- 
put field strengths, we describe in detail such a stability 
domain map in Ref. [26j], wherein the regions of stable 
switching, unstable regions that exhibit self-pulsing and 
chaotic dynamics are clearly identified. The bistable be- 
havior can be obtained for different combinations of the 
input /output fields while either keeping the other input 
field constant (i.e. OB response of \E\ n \ versus |-Ef"'| 
as well as l-Ef"*!; while holding \E™\ constant and other 
such combinations) or, varying the other input field with 
a specific functional dependence. 

The output levels, the threshold of switching, the range 
of input fields exhibiting bistability are all dependent on 
the parameters like the decay rates (both atomic and 
the cavity), the detunings (both atomic and the cavity) 
and the cooperative parameters corresponding to the two 
transitions. The novel effects arise due to the inclusion 
of feedback for both the fields that couple the three-level 
atoms. The two fields being complex (on interaction) 
their relative phases and strengths become critical factors 
in determining the dynamics. In the context of nonlinear 
dynamics the corresponding phase space, as well as, the 
parameter space is exceedingly large owing to fourteen 
different physical parameters (71,2, Ai, 2 , #1,2, K i,2, Ci,2, 
E{ n 2 ) all of which can be varied independently. Unlike 
in the OB with a control field (without feedback, as in 
Ref. [13]) it is the simultaneous interaction of the atom 
with both the fields and their independent feedback from 
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the cavities that self-consistently determine the output 
fields (both their phase and amplitude). Hence, with re- 
gard to the numerical implementation we express a note 
of caution as this setup does not permit an apriori choice 
of both the (complex) output fields. In the conventional 
computation of OB (in steady state and in the mean field 
limit) involving feedback for one field, one usually speci- 
fies the output field (which could be chosen to be a real 
value) and calculate the requisite unique (complex) in- 
put field. Such strategy is conventionally adopted due to 
the multivalued nature of the output field, because of the 
input-output relationship in form of the S~ shaped OB 
curve; however a given output field uniquely determines 
the input field. A generalization of a such a strategy 
fails because it is impossible to choose the amplitude and 
phase of both the output fields apriori. The input fields, 
the cavity fields and their interaction with the medium 
sclf-consistently determine the amplitudes and phases of 
output fields. Without taking into account the complex 
nature of the fields within the cavity one would miss out 
the hump like feature (discussed below) at low light lev- 
els [27j which is crucial and eventually transforms into 
the negative hysteresis as shown in Fig. [^Ja). In or- 
der to deal with such numerical constraints we used the 
Newton- Raphson method to obtain the steady state solu- 
tion of the nonlinear atom-field equations along with the 
boundary conditions. As expected an arbitrary choice of 
relative strength and phases of the output field does not 
necessarily correspond to a physically viable input field 
variation. 
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FIG. 2: (Color online) (a) The OB response for the U12 field 
for various values of the input field at Wi . (b) Enhancement of 
the population inversion on the |1) -H- \2) transition with (red) 
and without (blue) the feedback for the field at uj\ . Similarly, 
(c) and (d) indicate the absorption Im{p23,} and Im{pi2}, 
corresponding parameter values j/i = 25 and Ci = 1000 
(red curve), for blue dashed curve \yi\ = 0.5 and Ci = 0. 
The other relevant parameters are C\ = 1000, C2 = 100, «j = 
l, 7i = l,A i = O,0 i = O fori = 1,2 



In order to mimic a typical experimental situation we 
use real values for the input fields and compute the re- 
sultant complex output fields self-consistently along the 
different branches of the S-shaped OB response. We also 
note that in order to focus our attention to the new low 
input light OB regime in this paper, we have avoided 
the regimes involving multi-stability which is easily ob- 



tained in this double cavity OB system. There are two 
important features that we highlight here. The first one 
is the enhancement of the output field at L02 coupling the 
1 2} <H» 1 3) transition arising due to the creation of inver- 
sion in the |1) <->■ |2) transition which leads to the hump 
like response and eventually to negative hysteresis OB. 
The second feature is the novel switching characteristics 
wherein the output fields exhibit self-pulsing, moreover 
the two fields completely mimic each other in their tem- 
poral response. The hump like feature in the OB response 
indicates an enhancement of the field at LO2 arising due 
to suppressed absorption along the |2) «-» |3) transition. 
This occurs in the lower cooperative branch at low inten- 
sities of the ui2 input field, where initially the population 
is dominantly in the ground state |3). A large coopera- 
tive parameter C\ along the upper transition |1) o |2) 
leads to enhanced interaction with the u\ field resulting 
in the extraction of a significant fraction of the popula- 
tion into the excited states. Furthermore, as the upper 
state (|1)) population builds up there is a lowering of the 
influence of the field at u)2 on these atoms. This leads 
to the creation of inversion along the upper transition 
|1) «-» |2) i.e. pn > p22- This dynamics occurs due 
to the asymmetric choice of the cooperative parameters 
C\ > C2 at the cavity resonant condition Q\ = 9 2 = 0. 
However, similar dynamics can be obtained for compara- 
ble values of C\ and C2 in the bad cavity limit, i.e. with 
finite cavity detuning 62 The decrease in absorp- 

tion (/m{/523j) of the 102 field, accompanied by enhanced 
absorption (7ro{pi2}) at the lo\ field with and without 
the feedback for the uii field clearly demonstrates the re- 
liance of the enhancement of the field at L02 on cavity 
assisted inversion [see Figs. EIb),(c) and (d)]. With in- 
creasing incident field y\ , the above effects are enhanced 
and the hump like feature becomes quite exaggerated and 
eventually results in OB with a negative-hysteresis loop 
as indicated in Fig. [2a). A similar scenario is observed 
with increasing cooperative parameter C\ of the upper 
transition |1) o \2). 

We describe in detail the negative hysteresis loop aris- 
ing in this system (Fig. [3]). As described above the field 
at u>2 is enhanced, however, on further increase in the in- 
put field 2/2 the population from the upper states is drawn 
back into the lower levels |2) and |3) ultimately leading to 
large absorption of the o>2 field, and the output switches 
to the off- state (indicated by Ni in Fig. [3fa)). In the re- 
verse direction as the input intensity 2/2 is decreased the 
output field switches from an off- state to an on- state 
along N2 (different from N\ ) due to the large u>\ field al- 
ready circulating in the other cavity, thus encompassing 
within it a negative-hysteresis loop. This hysteresis is ex- 
actly opposite to the usual bistability (we denote as the 
positive hysteresis) wherein low input intensities leads to 
low output intensity and only for larger input intensi- 
ties the transition saturates (for zero atomic detuning) 
leading to the large output field and the corresponding 
reverse loop that encloses a hysteresis (indicated as P\ 
and P2 in Fig. |3fa)). This conventional positive hystere- 
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sis OB occurs at higher input field intensities where larger 
input field results in large output field and vice- versa. 




FIG. 3: (Color online) (a) The OB response for the u>2 
field which exhibits negative as well positive hysteresis (for 
Ci = 5000, C 2 = 20, \yi\ = 50). The color red (black) in- 
dicate stable(unstable) steady state response, and the arrows 
indicate the associated switching, (b) and (c) Critical slowing 
down associated with the negative hysteresis along Ni and 
./V2, respectively. The time evolution of fields are indicated 
by dashed lines for the operating points \y2\0p = 6.5(iVi) 
and 3.5(Af2) and the solid lines indicate operating points 
\y2\0p = 6.2(iVi) and 3.6(N 2 ). The corresponding threshold 
points are \y2\th = 6.13(iVi) and 4.05(A?2). The other relevant 
parameters are same as in Fig. [2] 

The nature of threshold points associated with the 
negative hysteresis are similar to the positive hystere- 
sis of the conventional OB and displays critical slowing 
down. The study of the time-dependent switching in- 
dicates that, as one chooses the operating point \y2\0p 
closer to the threshold of switching one observes an in- 
crease in the time required to switch to the steady state 
[Fig. |3fb-c)]. Note that this behavior occurs at all the 



four threshold points \y2\th associated with the switch- 
ing transitions N12 an d -Pi, 2- The variation of any other 
parameter results in shifting of the threshold point itself 
and thus the associated change in the switching times 
due to critical slowing down. 

Apart from such multicolored stable switching there 
is a wealth of dynamics this system can display. It ex- 
hibits periodic self-pulsing in the steady state for both the 
fields. A constant input intensity at both the frequencies 
wi and ui2 results in periodic output that largely mimics 
each other, and the periodicity can be controlled using 
the cavity parameters (such as Kx,2 and ^1,2)- These de- 
tailed nonlinear dynamical aspects are discussed in the 
companion paper [26) . Before we conclude, we would 
like to point out that there is nothing special about the 
two-photon resonant condition that we have used in the 
calculations presented here, similar results are obtained 
even in the non-resonant cases. The nonlinear dynamics 
can be obtained in a robust manner in a wide variety of 
regimes with appropriately chosen parameters. 

In conclusion, we have demonstrated a simple all op- 
tical double cavity OB system that exhibits negative as 
well as positive hysteresis. We self-consistently determine 
the amplitude and phases of both the output fields. A 
novel region of response involving a hump like feature in 
the S-shaped OB curves, as well as negative hysteresis 
is demonstrated at low input light levels. These effects 
are a consequence of the cavity induced inversion aris- 
ing from the simultaneous cooperative coupling at two 
different frequencies. The system also exhibits a range 
of nonlinear dynamical features such as self-pulsing and 
chaos, again at low input light levels. 
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